Parity non-conservation in thallium 



M. G. Kozlovj] and S. G. Porsev 
Petersburg Nuclear Physics Institute, Gatchma, 188300, Russia 



o 
o 

>>: 



Oh: 



o 

Oh: 



> 
o 
o\ 
o 
in 
o 



CZ3 



X 



W. R. Johnson 

Department of Physics, Notre Dame University, Notre Dame, IN 46556, USA 

(Dated: February 2, 2008) 

We report a new calculation of the parity non-conserving El amplitude for the 6pi/2 
6p3/2 transition in ^^^Tl. Our result for the reduced matrix element is SIpnc = —(66.7 ± 
1.7) ■ ilO~^^ {—Qw /N) a.u.. Comparison with the experiment of Vetter et al. [Phys. Rev. 
Letts. 74, 2658 (1995)] gives the following resuh for the weak charge of ^"^Tl: Qw/Q^ = 
0.97(±0.01)expt (±0.03)theor, where = -116.7 ± 0.1 is the standard model prediction. This 

result confirms an earlier conclusion based on the analysis of a Cs experiment that atomic PNC 
experiments are in agreement with the standard model. 

PACS numbers: PACS. 32.80Ys, ll.30.Er, 31.30.Jv 



I. INTRODUCTION 

The thallium atom is the second simplest atom where 
parity non-conservation (PNC) has been observed [Q, |j. 
The simplest and best understood such atom is cesium, 
where theory has the accuracy of 1% ||, ^ (see also more 
recent calculations [||, §,0]). That, together with the best 
experimental result for Cs ||] allowed one to obtain the 
most accurate result for the weak charge of the cesium 
nucleus li, ItII: 



0^(i33Cs) = 

QSM(133Cs) = 



-72.5 (3)cxpt(7)thcor, (1) 

-73.09(3), (2) 



where Eq. (|^) is the standard model prediction given in 
Ref. [^. Note, that these two values are in agreement 
with each other. Recently there has been some discussion 
of the accuracy of the theory for Cs; in Ref. [|oj the 
accuracy of the theoretical amplitude in Cs was estimated 
to be 0.4%, but at present there seems to be a consensus 
that it is closer to 1%. 

In this paper, we report a new calculation of the PNC 
i?l-amplitude for the 6pi/2 6^3/2 transition in thal- 
lium. Combining this new calculation with the most ac- 
curate experiment [Q , leads to a value of the weak charge 
of the 2°^T1 nucleus that can also be compared with the 
standard model prediction 

QwC°'T^l) = -113 (l)cxpt(3)thcor, (3) 

g^(205Tl) = -116.7(1). (4) 

Again, these two values are in agreement with each other. 
The atomic theory involved in both calculations is sim- 
ilar, as are the possible sources of theoretical uncer- 
tainty. That can possibly explain why the central values 
in Eq. (|^) and Eq. (|^) are shifted in the same direction 



from the standard model predictions (note that in both 
cases the theoretical uncertainty dominates the experi- 
mental one). 

The structure of this paper is as follows: In Sec. || 
we briefly describe the method of the calculation. More 
details can be found in Refs. |l|, ^ |l|. Results of our 
calc ulat ion and a discussion of its accuracy are given in 
Sec. III. Our conclusions are summarized in Sec. |v| . 



'Electronic address: mgk@MF1309.spb.cdu 



II. CALCULATIONAL DETAILS 

One important difference between thallium and cesium 
is that the former atom has three valence electrons above 
a relatively compact and rigid core, while the latter has 
only one. Correlations between three valence electrons 
are not small and can not be treated accurately by many- 
body perturbation theory (MBPT). In Ref. |jll|, it was 
suggested that MBPT be combined with the configura- 
tion interaction (CI) method for such atoms. The latter 
method is well suited to account for correlations between 
a few valence electrons, while the former method allows 
one to treat core- valence and core-core correlations. This 
combined CI-t-MBPT method was first tested in calcula- 
tions of energy spectra of Tl ||ll|, then for Ca, Sr, Ba, 
and Yb |Q . Later it was extended to the calculations of 
observable such as hyperfine structure constants [l^ ] and 
polarizabilities [^ . 

There is a significant difference in MBPT for one- 
electron atoms, such as Cs, and three-electron atoms 
such as Tl. For one-electron atoms one has to calculate 
only MBPT diagrams with one external line, while for 
many-electron atoms there are also diagrams with two, 
three, and more external lines. For combinatorial rea- 
sons the number of such diagrams grows rapidly with the 
number of external lines, making calculations for many- 
electron atoms much more complicated. Fortunately, the 
three-particle diagrams appear to be small for Tl 
If the three-particle diagrams are neglected, the eS'ective 
Hamiltonian for valence electrons is a two-particle oper- 
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ator, which not only drastically reduces the number of 
diagrams at the MBPT stage of the calculation, but is 
also essential for the CI stage of the calculation because 
the Hamiltonian matrix remains sparse. 

It was recently shown for Cs |^] that the Breit 

interaction correction is larger than earlier estimated . 
For this reason we have included the dominant magnetic 
part of the Breit interaction (the so-called Gaunt in- 
teraction) at all stages of the present calculation. The 
Dirac-Fock equations were solved in the Coulomb-Gaunt 
approximation for the ground state of the Tl"*" ion (the 
V^~^ approximation) with the help of the code described 
in Ref. @. 

We used different basis sets for the CI and MBPT parts 
of the calculation. That enabled us to improve the CI 
convergence without having an enormous Hamiltonian 
matrix. The core and valence orbitals were taken from 
the numerical solution of the Dirac-Fock equations on a 
radial grid. The CI basis set included virtual orbitals 
formed using the method suggested in Ref. jl^, while 
virtual states for the MBPT basis set were formed from 
_B-splines jl^. The Dirac-Fock Coulomb-Gaunt Hamil- 
tonian was diagonalized on both sets. Two different vari- 
ants of the MBPT basis set were used for calculating di- 
agrams for the effective Hamiltonian and for solving the 
RPA equations. The former included partial waves with 
1 = 5, while the latter included more orbitals but was 
restricted to / < 2. 

We use the Brillouin-Wigncr variant of MBPT. In this 
formalism, the effective Hamiltonian for the valence elec- 
trons is symmetric but energy-dependent. It was shown 
earlier |18|| that the accuracy of the theory can be im- 
proved by calculating the Hamiltonian at shifted ener- 
gies. The optimal valence energy for Tl was found to be 
-1.64 a.u.. 

For some observables such as polarizabilities, Stark- 
induced amplitudes and PNC amplitudes, one needs to 
sum over intermediate states or solve the corresponding 
inhomogeneous equations (Sternheimer ||l9| or Dalgarno- 
Lewis |20(] method). Here we apply the Sternheimer- 
Dalgarno-Lewis method to the valence part of the prob- 
lem as described in fO, Ell. 



III. RESULTS AND DISCUSSION 



A. Test Calculations of the Spectrum and 
Observables. 



To control the accuracy of our calculation of the PNC 
amplitude, we calculated as many different observables 
as possible. These results are presented in Tables | - 
[V. There are some differences with earlier calculations 
^n] , p^ . Firstly, we use the Coulomb-Gaunt potential 
here instead of the Coulomb potential used in those ear- 
lier calculations. In Table ||, we compare our theoretical 
Coulomb-Gaunt energies with the measured spectrum. 
We also give Coulomb energies for comparison. One can 



TABLE I: Spectrum of Tl for Coulomb and Coulomb-Gaunt 
potentials. The ground state three-electron binding energy 
l-Evail corresponds to the sum of the first three ionization po- 
tentials. For other levels, we give transition frequencies A. 
The eff'ective Hamiltonian for valence electrons is calculated 
for the energy iSvai = —1.64 a.u. 





Cl+MBPT 
C C-G 


Experiment 


6pi/2 -Eval (a.u.) 


2.0742 2.0720 


2.0722 |22j 



6p3/2 
7Sl/2 

7pi/2 

7p3/2 
6^3/2 



A (1/cm) 



6di 



'5/2 



7925 


7836 


7793 


26583 


26455 


26478 


34193 


34087 


34160 


35215 


35098 


35161 


36363 


36208 


36118 


36469 


36321 


36200 



TABLE II: Hyperfine constants for ^'^'^Tl (MHz). First row 
gives Dirac-Fock values and the following rows give various 
corrections described in the text. 







^6P3/2 ^"'"1/2 




^7p3/2 


^6^3/2 


^6d5/2 


DF 


17339 


1291 


7579 


1940 


187 


21 


9 


CI 


924 


-1369 


3799 


-102 


112 


-185 


391 




3428 


-45 


765 


331 


-56 


114 


-226 


Arpa 


959 


359 


1031 


103 


73 


5 


15 


A, 


-1071 


-31 


-269 


-92 


-9 


3 


-5 


Asbt 


-1389 


-161 


-75 


-113 


-19 


-19 


-8 


Atp 


1731 


120 


-22 


133 


4 


21 


7 


AsR 


209 


88 


-29 


14 


6 


-1 





Norm. 


-467 


-4 


-113 


-20 


-3 








Total 


21663 


248 


12666 


2193 


295 


-41 


183 


Theor." 


21760 


-1919 


12470 


2070 


195 






Theor.*" 


21300 


339 


12760 










Theor.'= 


21623 


264 


12307 


2157 


315 


-35 


184 


Expt. 


21311 


265 


12297 


2155 


309 


-43 


229 




c 


Ref.|^ 


Ref. 


|5| ^ Ref.y 







see that the Gaunt correction to the fine structure is 
rather large, especially for the ground state. Both split- 
tings 6pi/2-6p3/2 and 1pi/2~'^P3,/2 are significantly im- 
proved when the Gaunt interaction is included. The over- 
all agreement with experiment is also improved, though 
not as dramatically. 

Tables || and |H present results of our calculations of 
hyperfine structure constants and _Bl-amplitudes. In Ta- 
ble H, we list all corrections to the initial Dirac-Fock val- 
ues of hyperfine constants. The CI correction accounts 
for mixing of configurations. The ifeff correction cor- 
responds to the solution of the CI equations with the 
effective Hamiltonian and accounts for core polarization 
effects. The remaining corrections correspond to the dif- 
ferent terms in the effective operator of the hyperfine in- 
teraction |l^. In this calculation we included the struc- 
tural radiation {Asr) correction [ p6[ , which was omitted 
in . The most important contributions are associated 
with the RPA and the Brueckner {A„) corrections to the 
effective operator. Two other large contributions from 
subtraction (Agbt) and two-particle corrections (Atp) al- 
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TABLE III: Reduced matrix elements for iSl-amplitudes in 
L-gauge (a.u.). 









DF 


CI 


Total 


Thoor. 


Expt. 


6pi/2 




7Sl/2 


2.049 


1.863 


1.77 


1.72" 1.81(2)" 














1.78'' 








6d3/2 


2.722 


2.454 


2.30 


2.39" 


2.30(9)" 


6P3/2 




7Si/2 


3.966 


3.466 


3.35 


3.18" 
3.31'' 


3.28(4)" 






6^3/2 


1.633 


1.472 


1.40 


1.39" 


1.38(7)" 






6^5/2 


4.840 


4.292 


4.08 




4.0(2)" 














3.8(2)"* 


7Pl/2 




7Sl/2 
6^3/2 


6.618 
11.980 


6.152 
10.874 


5.96 
10.86 




5.87(8)" 


7P3/2 




7Si/2 


8.794 


8.252 


7.98 




7.88(11)" 






6^3/2 


5.395 


4.887 


4.90 










6(^5/2 


16.300 


14.799 


14.88 







" Ref.y '' Ref.|5[ " Refs.|7|, ^ Ref.|29| " Rcf.| 



TABLE IV: Polarizabilities of the 6pj-levels of Tl in a.u.. 
Atomic polarizability includes valence and core contributions. 
The column SCore accounts for the change of the core polar- 
ization due to the fact that some valence orbitals are occupied. 



ao(6pi/2) 

«o(6p3/2) 
Q2(6p3/2) 



Valence Core ^Core Total Expt. [pl[ 



43.47 
73.79 

-25.04 



6.23 
6.23 




-0.51 
-0.48 
0.06 



49.2 
79.6 
-25.0 



-24.2(3) 



most cancel each other. 

One can see that the Dirac-Fock approximation for the 
hyperfine constants works reasonably well only for 6pi/2 
and 7pi/2 levels. For some of the other levels, even the 
sign and order of magnitude of the constants are wrong. 
That affects the final accuracy of the theory, which is 
better than 2% for the npi/2 levels, about 3% for the 7s 
level, and worse for other small constants. It is actually 
quite surprising, that the CI-I-MBPT method gives rather 
accurate values even when the Dirac-Fock approximation 
fails completely. 

The MBPT corrections to the El amplitudes in the 
length gauge are much smaller than for the hyperfine 
constants and the final values are closer to the initial 
Dirac-Fock ones. Usually, the V-iorm of the El ampli- 
tude is less accurate for the optical transitions in neutral 
atoms. Therefore, we used the L-form for the calcula- 

Accurate experimental re- 



tions presented in Table III 



suits are available only for the npj —> Is amplitudes. For 
these four amplitudes the difference between theory and 
experiment is within 2%. Taking into account the small- 
ness of the MBPT corrections to the El amplitudes, we 
assume this to be the characteristic accuracy of the the- 
ory for El amplitudes. 

As an additional test of the theory, we calculated the 
tensor polarizability of the 6^3/2 level, which was accu- 
rately measured by Gould |Q (see Table IV). The con- 
tribution to the polarizability was calculated by solving 
the inhomogeneous equation for the effective operators in 



TABLE V: SlpNC amplitude for 6pi/2 - 
2°^T1 in the units i ■ 10"" (-Q^/iV) a.u. 



D3/2 



transition in 



right hand side: Hp\Gp^/2) Hp\&pz/2) 

CI -34.20 -29.88 

i^cft & RPA -7.26 -hO.Ol 

Brueckner -1-1.29 -1-1.12 

Subtraction -1-1.03 -1-0.77 

Two-particle -0.29 -0.53 

Struc. Radiation -0.04 -0.02 

Core sum -0.16 -1-0.06 

Subtotal -68.1 

Normalization -1-1.4 

Total -66.7 



the valence space. The core contribution was calculated 
as a direct sum of the RPA amplitudes over Dirac-Fock 
virtual orbitals. When calculating the polarizability of 
the core of the neutral atom, we should exclude occupied 
states from the sum over the intermediate states. The 
corresponding correction is given in the column JCore of 
Table The 3.5% difference between theory and ex- 
periment for 012(6^3/2) is in agreement with our estimate 
of the theoretical accuracy of the El amplitudes. 

We have omitted all three particle diagrams in the ef- 
fective Hamiltonian. These diagrams are strongly sup- 
pressed by the small overlap between valence and core 
orbitals [jll] . We calculated these diagrams for the small 
number of leading configurations and determined correc- 
tions to the energies and to the hyperfine constants. The 
valence energies changed by ^ 1Q~^ a.u. and corrections 
to the hyperfine constants were on the order of a few 
MHz, confirming that three particle corrections for Tl are 
well below the accuracy of the present calculation. We, 
therefore, neglected three-particle diagrams in all other 
calculations. 



B. PNC Amplitude. 

The PNC interaction of the electron with the weak 
charge of the nucleus has the form: 



Hp = 



9i 

2V2 



(5) 



where Gp = 2.2225 • 10"^^ a.u. is the Fermi constant of 
the weak interaction, 75 is the Dirac matrix, and p(r) 
is the neutron density of the nucleus. Calculation of the 
PNC amplitude is similar to calculation of polarizabili- 
ties. Firstly, we must determine the effective operator 
-ffp,cff for the PNC interaction of the valence electrons 
with the nucleus in the same way that we determined 
effective operators for the hyperfine interaction and the 
electric dipole moment operators. We must then solve 
the RPA equations for the operator (||), calculate the 
Brueckner Hp^^, subtraction TJpsbti two-particle Hp^p, 
and structural radiation -ffp^sR corrections. 
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Once the effective operators are formed, we solve two 
inhomogeneous equations: 



{P) 

6P3/2,™ 



- -ffp,off|^'6p3/2,m)- (7) 



Afterward, the PNC amplitude iJlpNc for the &pi/2 
6p3/2 transition can be readily calculated with the help of 
the (effective) electric dipole moment operator D = — er: 



i?lpNC = (*6p3/2ll^eff||«'^p,^,,) 
+ (*^pi/2ll^offll*6pi/2> 



(8) 



Alternatively, instead of (^) and (0) we can solve in- 
homogeneous equations with the operator ZJcff.o on the 
right hand side. The reduced matrix element of the PNC 
amplitude is then given by: 



El 



PNC 



= (-1)* 



3 1 1 

— m m 



(9) 



/2l-f^-P,cff|*6pi/2 



/2' 



Equations (||) and (g) should give identical results; a com- 
parison of results obtained using the two equivalent meth- 
ods can be used to check numerical accuracy of the cal- 
culations. On the other hand, it is easier to calculate the 
two-particle part of the operator Hp ^s using Eq. (|p) and 
the two-particle part of the operator Z?off using Eq. (@). 
The results of the calculation based on Eqs. (g) — (S) 
are given m Table 0. The inhomo geneous equations (n) 
and (|^) do not account for the sum over intermediate core 
states. This sum was calculated explicitly in the RPA ap- 
proximation and the corresponding results are also given 
in Table 0. 

a. Accuracy of the calculation. The principal the- 
oretical uncertainty in calculations for atoms such as 
Tl is associated with higher orders in the residual two- 
electron interaction. The effective operator method ac- 
counts for some important higher-order MBPT correc- 
tions. The valence-valence correlations are considered 
non-perturbatively within the CI method. In addition, 
by diagonalizing the effective Hamiltonian, we account 
for Brueckner and screening correction to all orders. The 
RPA equations for the PNC and the electric-dipole inter- 
actions effectively sum infinite chains of diagrams. How- 
ever, already in third-order MBPT, there are diagrams 
that are neglected here. 

It is known that for neutral atoms, the residual two- 
electron interaction is not small and it is very difficult 
to estimate the accuracy of different approximate meth- 
ods. Therefore, we made a number of test calculations 
described above. We have seen that, even for hyper- 
fine structure constants where the MBPT corrections 
are huge, the accuracy of the theory for the large con- 
stants is still about 2 - 3%, while for the El amplitudes 
the accuracy is better than 2%. According to Table 



the MBPT corrections to the PNC amplitude are rather 
small. Therefore, we estimate the error caused by the ne- 
glect of the higher-order terms in the residual interaction 
to be about 2%. 

The other possible sources of errors are QED correc- 
tions and the nuclear size effects. Below, we discuss both 
of them in some detail. 

b. QED corrections to PNC amplitude. The most 
important radiative corrections to the PNC interaction 
given in Eq. (|^) correspond to very small distances 
and should be calculated within the electroweak theory. 
These corrections lead to scaling of the weak charge of 
the nucleus. The leading-order (tree-level) value of the 
weak charge for ^os^pj jg. 

Qw{0) = -N + Z{l-4sm^0w)^ -117.9, (10) 

where N and Z are the numbers of neutrons and protons 
in the nucleus. The Weinberg angle 6w here is taken 
at the energy of Z-pole. Radiative corrections change 
Qw{0) by 1%; the resulting value, taken from Ref. Q, is 
given in Eq. (||). 

In addition to these radiative corrections, there are 
the QED radiative corrections specific to heavy atoms 
p2| , |3^ . These corrections are dominated by vacuum- 
polarization effects and for hydrogen-like ions with Z w 
80 were calculated to be 0.9% [|3|. For the valence 
electron of a neutral atom, screening of the vacuum- 
polarization potential by relaxation of the core should 
be taken into account. It was recently shown that for 
the Breit interaction this screening significantly reduces 
corrections to the valence amplitudes |^, Thus, 
the actual radiative corrections to the PNC amplitude in 
Tl could be few times smaller. 

Finally, there is the Breit correction to the PNC am- 
plitude. The magnetic part of the Breit interaction (the 
Gaunt interaction) is known to be larger than the ne- 
glected retardation part of the Breit operator [Q by a 
factor of three or more. As we pointed out above, the 
Gaunt interaction should be included self-consistently at 
all stages of calculations. It is difficult to isolate the 
Gaunt correction in the self-consistent approach. To do 
that, one must repeat all calculations in the Coulomb ap- 
proximation and take a difference. The calculations for 
Tl are rather time consuming, so for the Coulomb case 
we did only the CI calculation. At that level of approxi- 
mation, the Gaunt correction to the PNC amplitude was 
about —0.5%. We think that the overall Gaunt correc- 
tion may be up to two times larger. The retardation cor- 
rection, which is neglected here, should be at least few 
times smaller than the Gaunt correction and the total 
uncertainty from QED effects should be within 1%. 

c. Nuclear size effects. The PNC amplitude is sensi- 
tive to both proton and neutron distributions: the former 
determines electronic wave function inside the nucleus 
and the latter determines the weak charge distribution in 
Eq. (^). In our calculations we approximate the nucleus 
205T1 by a uniform ball of the radius 0.1334 • IQ-^ a.u. 
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This value corresponds to a root-mean-square charge ra- 
dius tn = 5.470(5) fm |3|. 

This model assumes that (i) proton and neutron radii 

of the nucleus are the same: r^"* = = tn, and (ii) 
the nucleus has a sharp edge. Generally speaking, both 
assumptions are incorrect. Therefore, it is important to 
estimate the errors associated with each of them. That 
was recently done in Refs. p6|, p7|: 



SEl 



PNC 



V 



-0.39 



„(«) JP) 



0.14 



N 



'N 



21 



25 H )2 



- 1. 



0.0377,(11) 
(12) 



The first two terms in (|llj) give the dependence of the 
PNC amplitude on the nuclear radius tn and on the dif- 



ference r' 



(») Jp) 



At present there are no accurate exper- 
imental data on neutron radii. However, theory predicts 



N 



'N 



0.1-0.3 fm 



The 



that for heavy nuclei rj^ 
parameter r] is defined so that it vanishes for a uniform 
distribution |3^, therefore, the last term in ( pT| ) describes 
the dependence of the PNC amplitude on the details of 
the nuclear distribution. For a real nucleus, i] is about 

0.1 0. 

Substituting into Eq. (pi]), we find that the correction 
to the PNC amplitude from the details of the nucleon 
distribution is 



SEl 



PNC 



El 



PNC 



-0.003(2). 



(13) 



This correction is seen to be small in comparison with 
other ones. 



These two calculations account for different correlation 
effects, but give very close results. In general, the am- 
plitudes for the allowed Ml-transitions are very stable 
and can be calculated quite reliably. The experimen- 
tal value of this amplitude follows from the measure- 
ments of the quadrupole amplitude E2 |39t] and the ratio 
X = uj/{2V3c)E2/Ml M: 



E2 = 13.29(3) 
X = 0.2387(40) 



Ml ^ 4.16(7) 10" 



(17) 



which is in agreement with theoretical results of Eq. ( p^ ) 

Using the theoretical value for the Ml amplitude ( p^ 
and the standard model value of the weak charge 
get the following result for the PNC rate TZ: 



we 



w 



-116.7) 



-15.2(4) • 10" 



(18) 



The experimentally measured PNC rate for the 6pi/2 
6p3/2 transitions is 



n 



-14.68(17), 
-15.68 (45), 



Vetter et al. 
Edwards et i 



(19) 



These experimental results formally contradict one an- 
other. Recently, Majumder and Tsai suggested that 
the discrepancy could be due to the different values of the 
parameter x used in the analysis by two groups. They ac- 



curately measured x ( 
the experimental resu 



ITD and used their value to rescale 
t from Rcf. [0 to find 



IV. CONCLUSION 



TZ = -14.71 (45), Majumder and Tsai [40]. (20) 



Our final value for the PNC amplitude of the 6pi/2 — > 
6p3/2 transition is in good agreement with the most ac- 
curate previous calculation of Dzuba et al. p6[ (t he ref- 
erences to earlier calculations can be found in [p8|): 



El 



PNC 



-(66.7 ±1.7) this work, 
-(66.1 ±2.0) Dzuba e< a/. 



ilO-ii(-Qw/iV) 
In the experiment the following ratio is measured 

ElpNC 



TZ = Ini 



Ml 



(14) 
(15) 



where Ml is the 6pi/2 ~^ 6733/2 magnetic-dipole transi- 
tion amplitude. We have calculated this amplitude in 
two different ways. Firstly, we did the CI+MBPT calcu- 
lation using the effective Hamiltonian. Secondly, we did 
a third-order MBPT calculation in which Tl was treated 
as a one-electron atom: 



Ml 



4.145- 10-^ a.u. 
4.149- 10-3 a.u. 



CI+MBPT-II, 
MBPT-III. 



(16) 



This scaled value is in agreement with the measurement 
||l|, where a nearly identical value of x was used, and all 
three values are in agreement with the theoretical result 
(|l|) for the Qw ^Qw^- We use the best experimental 
result |l| together with our calculation (^8|) to derive the 
experimental value of the weak charge for ^°^T1: 



Qw('°'Tl) = -113 (l)expt(3)thcor - (21) 
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